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We present a Hamiltonian derivation of a class of reduced plasma two-dimensional fluid models, 
an example being the Charney-Hasegawa-Mima equation. These models are obtained from the same 
parent Hamiltonian model, which consists of the ion momentum equation coupled to the continuity 
equation, by imposing dynamical constraints. It is shown that the Poisson bracket associated with 
these reduced models is the Dirac bracket obtained from the Poisson bracket of the parent model. 
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CN . I. INTRODUCTION 

Q, 

, Modelling plasma dynamics, with kinetic or fluid approaches, often amounts to investigating reduced, relatively 
tractable, models that capture the essential ingredients of the phenomenon under consideration while neglecting, for 
example, irrelevant spatial or temporal scales. Obtaining valuable reduced models provides a practical computational 
advantage when numerical simulations of a phenomenon have to be carried out. Reduced models are sometimes derived 

on : 

^vq from parent models through a well-defined reduction procedure that amounts to approximation directly at the level of 

voj ' 

the model equations, after having introduced some ordering based on physical arguments. Ideally, the reduced model 

t— I ■ 

should inherit some essential properties from the parent model. It has been recognized that the ideal part of such 

O ! 

parent models possess a Hamiltonian (although noncanonical) character, consisting of a Hamiltonian functional, which 

> : n 

k> ' can be identified as the total energyof the system, and a noncanonical Poisson bracket 111. Paradigmatic examples 

^ : nn n 

are the Vlasov-Maxwell equations [2|-|4| and the ideal MHD equations [5|, which both possess such a Hamiltonian 
structure. If one starts a derivation from a Hamiltonian parent model, then the final reduced model should also 
possess a Hamiltonian structure. If this were not the case, some faulty dissipation would enter the reduced model, 
and lead to qualitatively different interpretation of physical behavior. For instance, numerical simulations of a reduced 
model with spurious dissipative terms might converge to attracting states, which would be ruled out if the Hamiltonian 
structure were present. 

When the physical arguments invoked to reduce the parent model take the form of constraints on the dynamics, 



the method of Dirac brackets (see, e.g. Ref. 



6r|9(), provides a systematic method for obtaining a Hamiltonian reduced 



model from a Hamiltonian parent model. This method is general and works in the case of noncanonical Hamiltonian 
systems, which is the type of systems that arises when the models are formulated in terms of Eulcrian variables, 
fre quen tly used in plasma physics. Dirac brackets have proven to be useful for both finite and infinite-dimensional 



IOM12I] , for the derivation of balance models in geophysical fluid dynamics 



13, 



14j , as a numerical (simulated annealing) 

. The 



method for calculating vortex states |15j , and also for describing the dynamics of fluids with free boundaries 
purpose of the present paper is to show how a class of reduced models, which will be Hamiltonian by construction, 
can be obtained from a Hamiltonian fluid parent model. Relevant examples of reduced models belonging to this class 
include the Charney-Hasegawa-Mima equation 17 



18j . describing the propagation of drift waves in plasma, and the 
Euler equation for an incompressible fluid. 

Here a slab geometry with Cartesian coordinates is adopted and the dynamics of a plasma is confined to a plane 
transverse to an imposed magnetic field. Given a uniform and constant magnetic field B = Bz, the parent model 
describes the ion fluid dynamics in terms of its density and velocity fields. This model has a noncanonical Hamiltonian 
form 5|. The question addressed is how to derive reduced models describing the evolution of a few fields (for instance, 
the density or the electrostatic potential), given some physical constraints dictated by experimental relevance, while 
preserving the Hamiltonian character of the parent model. In order to perform the reduction, two constraints are 
considered: First an incompressibility assumption on the ion fluid, and second a relation between the density and 
the electrostatic potential (or equivalently of the streamfunction of the incompressible part of the ion velocity field). 
This last constraint can eventually be ascribed to a relation that couples the ion and the electron fluid through a 
quasi-neutrality assumption. We show using Dirac brackets that the Poisson bracket of the reduced fluid models can 
be constructed and the noncanonical Hamiltonian structure of these reduced models are recovered in a systematic 
way. 

The paper is organized as follows: In Sec. [IT] we recall the Hamiltonian structure of the parent model and introduce 
new dynamical variables that are particularly convenient for the problem under consideration. In Sec. lIIII the essential 
elements of Dirac 's theory of constrained Hamiltonian systems are briefly reviewed and the explicit derivation of the 
Dirac bracket for the reduced models is carried out. In Sec. lIVI we focus on the specific examples of the Euler equation 
for an incompressible fluid and the Charney-Hasegawa-Mima equation. We also discuss here the model derived by 
Terry and Horton [l9(. Conclusions are drawn in Sec. tVl 
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II. EXPANSION OF THE ION FLUID DYNAMICS AROUND EQUILIBRIUM 



We start the derivation from a parent model with two dynamical equations: one describing the transverse dynamics 
of the ion velocity field v(x, y, t) and the other describing the dynamics of the ion density field n(x, y, t): 



v + (v • V)v = -Vp + v x B, 
h = — V • (nv), 



(1) 
(2) 



where the dot indicates the partial derivative with respect to time t. We have used units such that the ion mass is 
M = 1, its charge e = 1, and the amplitude of the magnetic field B = 1. The total energy of the ions, given by the 
sum of their kinetic energy plus the potential energy provided by the electrostatic potential ip, is a conserved quantity 
that is also a good candidate for the Hamiltonian of the system of Eqs. ([IE]), viz. 



H\n,v] = J d 2 x 



The dynamics is determined by the Poisson bracket 5 



21] 



v 

i — 
2 



{F,G} 



dfx 
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- nip 
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n ) 





(3) 



(4) 



where we denote the functional derivatives of a given observable -F[n,v] by subscripts, i.e. F v = SF/Sv and F n — 
SF/5n. In our context we assume that the electrostatic potential ip is determined by the dynamics of the electrons 
which leads to a function <p(n e ), where n e is the electron density. From the quasi-neutrality condition, n — n e , the 
Hamiltonian becomes 



H[n,v] = J 



d 2 x 



ip(n) 



(5) 



where ip' '(n) = <p(n). We notice that for an external potential <p, we obtain ip(n) — nip as in Eq. ([3]). Another example 
is obtained by neglecting the inertia of the electrons so that their density obeys a Boltzmann law n e — no exp p, where 
no = no(x,y) = 1 — A(a;, y) is the electron density at equilibrium, given by a constant plus a space dependent part 
that we will later assume to be small. Hence, p(n) = log(7i/(l — A)) and consequently ip(n) = n(log(n/(l — A)) — 1). 
We perform the following change of variables (n, v) i->- (h, <fi, D) defined by 

n = n, 

Acj) = z • V x v, 
AD = V • v , 



where A denotes the Laplacian. For simplicity, we use n instead of h in what follows. In terms of the new variables 
(n,<p,D), Hamiltonian ([5]) becomes 

•|V0| 2 + |V£| 2 



H[n,(f>,D] = J d 2 x 



+ [</>, D] +V(n) 



(6) 
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where [f,g] — z • V/ x Vg, and the bracket (j4]) becomes 

{F,G} = J d 2 x(F n G D -F D G n + ^^{[A- x F^A- x G^} 

+ [A- 1 F D ,A- 1 G D ] + VA^Fd • VA- X G - VA" 1 ^ • VA^G^ 

We first assume that the variables evolve slowly with time, which is equivalent to adding a prefactor of 1/e to 
the Hamiltonian, and we introduce an e-ordering for the dynamical variables. The hypothesis is that the system of 
interest is near an equilibrium state whose spatial variations are of order e: 

n = 1 + enx, 
<f> = e<f>i, 
D = eD l . 

The expansion of the Poisson bracket up to order 0(e ) is given by 

{F,G} = i J d 2 x(F ni G Dl -F Dl G ni -F Dl A- 1 G^+F^ 1 A- 1 G Dl ) 

-i J d 2 x (F ( f >1 A~ 1 £A~ 1 G ( t >1 + F Dl A^ 1 CA^ 1 Gd 1 — F Dl A~ 1 AA~ 1 G c t> 1 + F ctll A~ 1 AA~ 1 GD 1 ) , (7) 

with the linear operators C and A defined by 

Cf=[A<f >1 -n 1 ,f], 

Af = -V • ((A^i - m)V/) , 

where we omit their dependence on Ac/>i — n\. Observe, C is anti-self-adjoint {C' = —£), while A is self-adjoint 
(At = A). 

In the next section, we impose constraints on the Poisson bracket ((TJ) and compute the associated Dirac bracket. 

III. DIRAC BRACKETS 



First we recall few basic facts about Dirac brackets in infinite dimensions. If we impose N Eulerian constraints 
3> Q (x) = for a = 1, . . . , JV on a Hamiltonian system with a Hamiltonian H and a Poisson bracket {-,-}, the Dirac 



bracket is obtained from the matrix C defined by the Poisson bracket between the constraints 

C ab (x,x') = {$ tt (x),$ 6 (x')}, 

where note C Q ^(x, x') = —Cp a (x',x). The Dirac bracket is defined by 

{F,G}, = {F, G} ~Jd 2 xJ d 2 ^{F,$ a (x)}C^ 1 (x,x'){^(x'),G} ) (8) 

where the C~i(x, x') are defined by 

/ ( iVC- /9 1 (x,x')^ 7 (x , ,x") = y d 2 x , C a/3 (x,x')C^ 7 1 (x',x") = <5 a7 <5(x - x"), 

which implies C~^(x,x') = -C^(x',x). 

Dirac obtained ([5J from a modified Hamiltonian with Lagrange multipliers associated with each constraint: 

H' = H + J d 2 xM Q (x)$ Q (x). 

The computation of the dynamical equation associated with this new Hamiltonian gives 

{F,H'} = {F,H} + J d 2 x» a {F,$ a } + J d 2 x$ a {F,[i a }, 

which is equal to 

{F,H'}n{F,H}* = {F,H} + J d 2 xfi a {F,^ a }, 

where the symbol w means equality after the constraints are imposed. The coefficients /i Q are obtained by demanding 
that $ Q are Casimir invariants of the bracket {•, •}», which leads to the definition of the bracket ©. 

We impose two constraints on the dynamics. The first one is incompressibility, which translates here into 

$ x (x) = £>! = {), 

and the second one is an assumption relating the electron density to the electrostatic field. This assumption takes the 
form 

$ 2 (x) = ni-JV(^i) = 0, 

where N is, in general a nonlinear pseudo-differential function of <f>\ , i.e. a function of 4>\ and its derivatives to arbitrary 
order. The functional derivative of $2 with respect to 4>\ is given by 

— — — = -NH x -x), 
<5<M X ) 



where N is the Frechet derivative of N defined by 



NSfa = -y-JV^i +e8<t>x) 
de 



The Poisson brackets between the constraints are given by: 



6=0 



Cn(x,x') = -iA- 1 £A- 1 J(x-x'), (9) 
C 12 (x, x ') = -i fl - A- 1 ^ + eA-^A- 1 ^) <5(x - x'), (10) 



C 21 (x', x) = i fl - 7VA- 1 + eiVA^AA- 1 ) <5(x - x'), (11) 
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C* 22 (x,x') = -iivA-^A^jV^fx-x'). (12) 
e 

We take the convention that the linear operators act on the first mentioned independent variable, e.g., 7V5(x — x') 
involves the derivatives with respect to x, while £<5(x' — x) involves the derivatives with respect to x'. The elements 
of C _1 are determined by the following four relations: The first two determine C-j^x, x') and C^^x, x'): 

J d 2 x> [C 11 (x ) x')C 1 - 1 1 (x',x") + C 12 (x,x')C 2 i 1 (x',x")] = 5(x-x"), 
J d 2 x' [C 21 (x,x')C n 1 (x',x") + C 22 (x,x')C 2 l 1 (x',x")] =0, 

which upon using Eqs. (IWT^j) become 

-eA^CA- 1 ^ 1 - (1 - A- 1 ^ + eA _1 AA _1 7V t )C 21 1 = e 2 5(x - x'), 
(1 - NA- 1 + dVA~ 1 AA -1 )C , ii x - eN A' 1 CA- 1 = 0. 

The expansion of the solution of the above equations, up to order 0(e 4 ), is given by 
Cfj 1 (x,x / ) = -e 3 (l - JVA-^-^A-^A-^Cl - A- 1 7V t )-M(x - x') + 0(e 4 ), 

O^fcx') = -e 2 (l - A- 1 iVt)- 1 5(x-xO + e 3 (l - A- 1 iVt)- 1 A- 1 AA- 1 A^t(l - A- 1 iVt)- 1 ( 5(x-x') + 0(e 4 ). 
In the same way we obtain the expansions for (x, x') and CjTj^x, x'): 

/" dV [C 22 (x,x / )CS- 2 1 (x / ,x // ) + C r 2 i(x,x')C , i- 2 :l (x / ,x")] = 5 ( x - 
/" dV [^(x^OC^^x'^'O + Ci^x^O^^x',^')] =0, 

which upon using Eqs. (|9|)-(fT2j) leads to 

C£(x,x!) = -e 3 (l - A" 1 iVt)- 1 A- 1 £A- 1 (l - ^A" 1 )" 1 .^ - x') + 0(e 4 ), 

Cf^x') = e 2 (l - NA' 1 )-^* - x') - e 3 (l - iVA-^^iVA^AA-^l - ^A -1 ) -1 ^ - x') + 0(e 4 ). 



Now, using the above and the following expressions: 



{F, $i(x)} = ? (F ni + A- X F 01 - eA^AA- 1 ^ + eA^A" 1 ^) , 

{F, $ 2 (x)} = -i ((1 - iVA- 1 )^ + eNA^CA^F^ + eNA^AA^Fo, 



we get 



NA'^^NA^CA^G 



<Pi i 



J d 2 x J d 2 x'{^,$ 1 (x)}Cr i 1 (x,x'){$i(x'),G} = 

i /" rf^^ +A" 1 F 01 )(l-iVA- 1 )- 1 iVA- 1 /:A- 1 iVt(l-A- 1 iVt)- 1 (G ni +A- 1 G 01 ), 
yd 2 x y ( i 2 x'{^$ 1 (x)}Cr 2 1 (x,x'){$2(x'),G} = l Jd 2 x(F ni +A- 1 F 4>1 )G Dl 

-i y d 2 x ((A^AA" 1 !^ - A-^A-^dJG^ - (F ni + A" 1 J F <Al )(l 
y d 2 * y d 2 x'{ J F,$ 2 (x)}G 21 1 (x,x'){$i(x'),G} = -^ J d 2 x F Dl (G ni + A^G^) 

+i y d 2 ^ 1 (A- 1 AA- 1 G 01 -A- 1 £A" 1 G i3 J-iVA- 1 £A- 1 ^ 1 (l-A- 1 iVt)- 1 (G ni +A- 1 G^), 
y d 2 a; y d 2 x'{ J F 1 ,$ 2 (x)}G 22 1 (x,x'){$ 2 (x'),G} = 1 y cPxF^^CA^Gd^ 

By summing these contributions and subtracting the result from the Poisson bracket (0, as per ©, we obtain the 
following Dirac bracket: 

{F, G}* =-- / ^(A-JV+r 1 ^ +7V t F„ 1 )£(A-7V t )- 1 (G 01 +iV+G ni ). 



If we define F for any functional ni, Di] as 



F[0r] =F[0 1 ,iV(^ 1 ),O], 



the functional derivative of F with respect to <j)\ is given by 



hi =F <t>1 +N*F ni . 



Therefore, the Dirac bracket becomes 



{F, G}» = ^ j d 2 x{A^-N(^)) (A-iVt)- 1 ^ 1 ,(A-7Vt)- 1 G 01 



(13) 



Finally, if we perform the following change of variables 



qi = Afa-N&i), 



the Dirac bracket achieves the compact form 

{F,G}* = - e J d 2 x qi [F qi ,G qi ], (14) 



as obtained in Ref. 21|. In particular, we notice that the Poisson bracket (fl4"l) does not depend explicitly on the 
constraint function TV. Therefore, all the models with constraints of the form m = N((p±) share the same Poisson 
bracket when expressed in terms of the generalized vorticity variable q\ . 

IV. HAMILTONIANS 

Having obtained a common Poisson bracket for our class of systems, we now obtain Hamiltonians for various cases. 

A. 2D Euler equation 

When the electrostatic potential is set to zero, the expansion of Hamiltonian ([5]) is given by 

Hi[ni,vi] = | j ' d 2 xv\, 

which in terms of the variables <j)\ and D\ is 

Hi = \J d 2 x (|V^i| 2 + |V£>i| 2 ) . 

As constraints we choose a constant density, i.e., N — 0, and incompressibility, i.e., Di = 0. The resulting Dirac 
bracket is given by 

{F, G}» = i J AA^fA-^.A- 1 ^]. 
This corresponds to the Hamiltonian and the Poisson bracket of the two-dimensional Euler equation, which was given 



23 



24| in terms of the vorticity variable qi = A<pi : 

H t = -- J d 2 xq 1 A- 1 q u 



and 



{F,G}« = J d 2 x gi [F qi ,G qi }. 
Since e does no longer play a role in the discussion, we drop it in what follows. 



B. Charney-Hasegawa-Mima equation 

Now we assume a linear adiabatic response of the electrons. Upon expanding ip(n) — log(n/(l — eA)), this gives 
N(<j>i) — (f>i — A and N = 1. The Dirac bracket becomes 

{F, G}* = J d 2 x (A0! - 4>i + \) [(A - l)- 1 ^ , (A - l)- 1 ^] , 



given by 



21] 



which is exactly the Poisson bracket found in Ref. 20| and derived in Ref. 2l| . The expansion of Hamiltonian ^ is 



Hi = \J d 2 x (|V0!| 2 + IV^I 2 + n 2 + 2Xm) , 

which becomes 

Hi = \J d 2 x (|V^| 2 +^ 2 ) , 

when introducing the constraints n\ = N(<f>i) and D\ = 0. 

We notice that for N{4>i) — 4>ihl where rj is small, we recover the Hamiltonian structure of the asymptotic model 
given in Ref. 22]. 



C. General setting 

We consider a Hamiltonian which is given by the sum of a kinetic energy plus a potential part 

W 2 



H = \ d l x 



^- + V[N(<H)], 



(15) 



where 



V[n] = / d 2 xip{n) , 



(16) 

is a quite general functional of the density n and tp is arbitrary. With this Hamiltonian and the bracket of (|13|) or 
(|14p , the equation of motion for <f>\ is given by 

(A - Njfa = - [A0! - N{<tn), (A - N^H^ 

where 
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This equation which originates from a Hamiltonian system is valid for any constraint N and any potential V. As 
stated above, all these models share the same Poisson bracket (in the generalized vorticity variable). They differ in 
their Hamiltonian. For V given by Eq. (fT())) and for %jj such that 

1 p'(N(4>t)) = <l>i, 

this equation simplifies, 

(A-JV)0i = [A0i-JV(0 1 ),0 1 ] 1 

which is verified for the two cases described above, the 2D Euler and the CHM equations. As a consequence of our 
ordering, this choice of ip is equivalent to the condition mentioned in SeclU namely that ip'(n) — <p(n). 

The question now becomes whether or not it is possible to construct a potential V for a dynamical equation or an 
electron response of interest. A rigorous answer is determined for the dynamical equation for the electric potential 



ip generated by the electrons. For example, in the nonlinear drift wave model derived by Terry and Horton 19( (see 



also 



25 



26j). the electric potential is not only a function of n but of all its derivatives, namely, 

ip(n) = N- 1 (n), 

where JV -1 is the inverse of the relation n = N(<p). In particular, N(ip) = (1 + O a )ip — A, where O a = <5o( c i + A)9y 
is an anti-self-adjoint operator. For such cases, the potential function ip with ip' = <p cannot be constructed. General 
potential functions of the form of that of (|15|) . such that 



V n = ip(n) 

can only be solved if is a self-adjoint operator. Thus the Terry-Horton model is not a Hamiltonian model in our 
class. 

V. CONCLUSIONS 

We have shown how the theory of constrained Hamiltonian systems, developed by Dirac, provides an effective way 
to include constraints into a Hamiltonian system while preserving the Hamiltonian character. In particular, we have 
applied this theory to a case relevant for plasma physics, by adopting as the parent model a fluid system consisting 
of the ion momentum equation and the continuity equation, where the electron dynamics determines the electrostatic 



11 

potential as a function of the density. We derived the Poisson bracket for the Euler and the Charney-Hasegawa-Mima 
equations as the Dirac bracket of this dynamics, obtained by considering an incomprcssibility condition and a general 
relation between the density and the velocity field as constraints. Given the bracket thus obtained, different choices 
for the Hamiltonian functional lead to different models, and in this way a general class of systems was derived, all of 
which preserve the constraints and are Hamiltonian by construction. 

We emphasize that the technique presented here is general, and can be applied to derive new reduced Hamiltonian 
models from more general parent models, models for which the Hamiltonian structure becomes directly available. 
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